Australian  Government 
Department  of  Defence 

Defence  Science  and 
Technology  Organisation 


Approximation  of  Bit  Error  Rates  in  Digital 

Communications 

Graham  V.  Weinberg  and  Sharon  Lee 

Electronic  Warfare  and  Radar  Division 
Defence  Science  and  Technology  Organisation 

DSTO-TN-0761 


ABSTRACT 


This  report  investigates  the  estimation  of  bit  error  rates  in  digital  communi¬ 
cations,  motivated  by  recent  work  in  [6].  In  the  latter,  bounds  are  used  to 
construct  estimates  for  bit  error  rates  in  the  case  of  differentially  coherent 
quadrature  phase-keying  with  Gray  coding  over  an  additive  white  Gaussian 
noise  channel.  By  analysing  Marcum’s  Q-Function,  which  is  an  integral  part 
of  bit  error  rate  expressions,  we  derive  more  direct  methods  of  estimation, 
including  least  squares  and  truncated  series  approximations.  Accurate  and 
efficient  estimates  for  bit  error  rates  are  then  proposed. 
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Approximation  of  Bit  Error  Rates  in  Digital 
Communications 


EXECUTIVE  SUMMARY 


DSTO  is  a  sponsoring  partner  in  the  Hypatia  Scholarship  scheme  for  mathematically 
talented  women,  and  as  such,  through  Electronic  Warfare  and  Radar  Division’s  Microwave 
Radar  Branch,  sponsored  the  second  author  of  this  report  to  participate  in  a  short  research 
project  through  the  Summer  Vacation  Scholarship  Program.  Named  after  the  famous 
female  mathematician  Hypatia,  the  scheme  provides  both  financial  and  mentoring  assets 
to  encourage  women  to  pursue  their  interests  in  the  mathematical  sciences.  The  work 
presented  here  is  a  report  on  this  project,  jointly  undertaken  by  Graham  V.  Weinberg  and 
Hypatia  Scholar  Sharon  Lee,  over  the  2006/2007  Summer  Period. 

This  project  involves  estimation  of  the  Marcum  Q-Function,  which  is  an  important  tool  in 
digital  signal  processing.  R  is  of  interest  to  both  the  radar  and  communications  research 
communities,  and  has  been  investigated  by  the  first  author  quite  extensively.  Here  we 
examine  bit  error  rate  estimation  in  digital  communications,  which  is  intimately  related 
to  this  function.  We  show  that  a  method  applied  in  a  recent  publication,  which  uses 
bounds  to  estimate  bit  error  rates,  can  be  improved  considerably  by  using  more  direct 
techniques  of  estimation. 

This  work  is  relevant  to  the  long  range  research  efforts  into  radar  detection  issues  associated 
with  Task  AIR  04/206,  EWRD  Support  for  AP-3C  E/LM2022  Radar  System.  Although 
focusing  on  a  communications  application,  the  results  transfer  directly  to  the  latter.  The 
technique  examined  here  will  be  useful  for  engineers  and  scientists  looking  for  efficient  and 
accurate  approximations  for  intractable  integrals. 
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1  Introduction 


This  report  examines  the  estimation  of  bit  error  rates  (BERs)  in  digital  communications 
[1].  Specifically,  we  will  investigate  the  recent  work  of  [6]  on  using  bounds  to  construct 
approximations  for  differential  quaternary  phase  shift  keying  (DQPSK)  transmission  with 
Gray  coding  over  an  additive  white  Gaussian  noise  channel  (AWGNG)  [14],  In  [6],  an 
estimate  of  this  BER  is  constructed  by  averaging  a  lower  and  upper  bound.  We  show 
that  more  direct  methods  can  be  applied  to  estimate  the  BER,  and  in  some  cases  more 
accurate  results  are  obtained. 

The  BER  is  a  fundamental  performance  measure  of  a  system,  quantifying  the  reliability 
or  integrity  of  a  received  signal  [1].  The  instantaneous  BER,  for  many  practical  com¬ 
munication  systems,  in  particular,  wireless  communications  systems,  can  be  written  as  a 
function  involving  the  standard  Marcum  Q-Function  ([11],  [12]  and  [13]).  This  famous 
function  has  received  much  attention  in  the  digital  signal  processing  literature,  due  to  its 
intractability.  Hence,  many  estimation  schemes  have  been  proposed,  employing  techniques 
such  as  Adaptive  Simpson  Quadrature  (ASQ)  [10],  Taylor  Series  approximations  [16],  the 
construction  of  lower  and  upper  bounds  ([4]  ~  [9]  and  [19])  and  the  Monte  Garlo  scheme 
([23]  and  [24]). 


1.1  Bit  Error  Rates  and  Marcum’s  Q-Function 


We  specialise  our  analysis  to  the  case  considered  in  [6],  as  described  above.  In  this  scenario, 
the  instantaneous  BER  is  described  by 

BER(7la,6)  :=  q{a^,b^)  -  hiab^),  (1) 

where  constants  a  and  b  depend  on  the  modulation/demodulation  format,  and  7  is  the 
total  instantaneous  signal  to  noise  ratio  per  bit,  and  Q(a,/3)  is  the  standard  Marcum 
Q-Function,  defined  by 


POO  1/2  2  \ 

Q{a,P)  :=  /  )lo{ax)dx,  (2) 

and  Io{-)  is  the  modified  Bessel  function  of  order  zero.  In  the  case  considered  in  [6], 

a  =  ^  and  b  =  \/2  +  ^72,  so  that  b  >  a.  The  key  to  estimating  (1)  is  to  construct 

good  approximations  of  (2). 

Using  lower  and  upper  bounds  on  (2),  derived  in  [5],  one  can  show  that 

BER(,i„, .) .  vf  {/l(»  - «)) .  (3) 

where  Erfc  is  the  complementary  error  function  (see  [6]  for  details). 
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Using  asymptotic  approximations,  [6]  also  proposed  the  following  approximate  expression, 
for  large  values  of  7: 


BER(7|a,6) 


b  +  a  - 
e 


VS-Kab  b  —  a 


(b-ar 


V^+1  1  ^[-(2-v/2h]_ 


(4) 


Although  (3)  and  (4)  produce  good  results,  we  will  investigate  a  more  direct  approach 
to  the  estimation  of  the  BER.  The  Least  Squares  (LS)  approach  is  an  interpolation  tech¬ 
nique,  which  may  be  applied  usefully  to  facilitate  the  estimation  of  (2).  Additionally, 
we  will  consider  Taylor  Series  approximations  applied  directly  to  (2).  Before  considering 
estimators,  we  introduce  a  number  of  useful  representations  of  the  Marcum  Q-Eunction. 


1.2  Integral  Representations  of  Marcum’s  Q-Function 


A  number  of  new  representations  of  (2)  have  been  derived  in  [22].  In  the  spirit  of  [4], 
these  convert  the  Marcum  integral  to  one  on  a  finite  domain,  with  penalty  terms  added. 
Converting  the  Marcum  integral  (2)  to  one  on  a  finite  domain  has  the  potential  to  improve 
the  estimation  process,  since  it  is  somewhat  easier  to  estimate  an  integral  on  a  finite 
domain. 


It  can  be  shown  that 


Q(a,6)  =  -  [l  +  e-“'/o(a2) 


+  xe  2(*^+“^)/Q(ax)dx. 


(5) 


Additionally,  by  an  application  of  the  symmetry  relation  [15]  of  the  Marcum  Q-Eunction 
to  (5),  one  can  derive 


Q(a,  b)  =  ^  [1  -  +  e-H“"+^")lo(a6) 


—  /  xe 


J  a 

Eurther  details  of  the  derivation  of  (5)  and  (6)  can  be  found  in  [22]. 


(6) 


It  is,  of  course,  not  difficult  to  reduce  (2)  to  an  integral  over  a  finite  domain.  Clearly  we 
can  write 


rb 

Q{a,b)  =  1  —  Io{ax)dx. 

Jo 


(7) 


Also,  one  can  apply  the  Marcum  symmetry  relation  [15]  to  (7)  to  produce 

Q(a,6)  =  e-K“"+^")/o(a6)-  /“ (8) 

Jo 
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We  will  construct  a  number  of  estimators  based  upon  (5),  (6)  and  (8),  using  the  Least 
Squares  Method. 


1.3  Least  Squares  Method 


The  Least  Squares  Method  (LSM)  [3]  fits  a  smooth  curve,  with  minimum  error,  to  a  given 
set  of  data  points.  Error,  in  this  case,  refers  to  the  sum  of  the  squares  of  the  errors 
(SSE)  or  the  residuals  of  the  points  from  the  curve.  To  fit  a  polynomial  of  degree  n, 
Pn{x)  =  On®"'  +  +  •  •  •  +  oix  +  uq,  to  a  Set  of  TTi  data  points,  where  m  >  n  +  1, 

the  coefficients  of  the  polynomial  are  determined  such  that  the  SSE  is  minimised,  where 


SSE  =  J2[Pnixk)-yk]\  (9) 

k=0 

and  {xk,  Uk),  k  =  1, . . .  ,m,  are  the  data  points. 


The  minimum  SSE  is  obtained  by  taking  the  partial  derivatives  of  (9)  and  equating  them 
to  zero.  The  resultant  set  of  equations  are  known  as  the  normal  equations  [3].  The 
coefficients  are  obtained  by  solving  the  normal  equations. 


m  m 

(y.\  ^  ^  Xji  “h  *  *  *  H“  Otn 
k=0  k=0 

mm  m 

Xk  +  cxi'^  xl  - h  ttn  X] 

k=0  k=0  k=0 


Y.yk 

k=Q 

m 

^XkVk 

k=0 


mm  m 

aoY.xl  +  aiY,  +  •  •  •  +  an  ^ 

k=0  k=0  k=0 


m 


^kVk- 

k=0 


(10) 


1.4  Contributions  of  this  Report 


This  report  explores  the  idea  of  applying  LS  approximations  to  evaluate  the  BER  Function 
(1),  achieved  by  estimating  the  Marcum  Q-Function  (2).  Instead  of  applying  bounds  to 
estimate  the  Marcum  Q-Function,  as  in  [6],  we  show  that  more  direct  methods  can  be 
applied  to  produce  estimators.  Specifically,  six  estimators  can  be  constructed  from  the 
finite  interval  representations  of  the  Marcum  Q-Function.  In  particular,  three  of  the 
estimators  use  an  exponential  function  with  a  quadratic  polynomial  power  to  estimate  the 
modified  Bessel  function.  The  second  group  of  three  fit  a  polynomial  of  degree  5  to  the 
entire  integrand.  These  six  estimators  are  compared  to  the  results  obtained  from  ASQ. 
Speed  and  accuracy  performance  are  also  analysed.  We  attempt  to  identify  an  optimal  LS 
estimator. 
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In  addition,  we  examine  a  Taylor  Series  approximation  of  the  BER  fnnction.  Its  perfor¬ 
mance  is  also  compared  with  ASQ  and  the  best  performing  LS  estimator. 
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2  Numerical  Approximation  Schemes 


A  number  of  estimators  of  the  BER,  based  on  the  LS  scheme  and  Taylor  Series  approach, 
are  now  introduced.  The  LS  estimators  are  constructed  by  substituting  functional  ap¬ 
proximations  to  the  integrands  in  (2),  (5),  (6)  and  (8).  In  some  cases,  the  modified  Bessel 
function  is  approximated  by  an  exponential  function,  while  in  others,  we  apply  a  polyno¬ 
mial  approximation  to  the  entire  integrand.  We  also  introduce  a  single  estimator  based 

upon  a  truncated  Taylor  Series.  Throughout  we  assume,  as  in  [6],  that  a  =  ^2  —  -v/2  and 

h  =  ^2  -|-  -v/2,  and  hence  b  >  a. 

The  key  problem  in  terms  of  estimation  of  the  Marcum  Q-Function  is  the  presence  of  the 
modified  Bessel  function  in  the  integrand.  Hence,  we  apply  functional  approximations  to 
eliminate  it  from  the  integrand,  to  facilitate  integration.  Many  of  the  proposed  functional 
approximations  suggested  below  have  been  produced  by  considering  suitable  fits  to  the 
Bessel  function  by  an  appropriate  polynomial.  This  will  be  done  by  using  Least  Squares 
fits. 


2.1  Least  Squares  on  a  Finite  Interval 


To  begin,  we  consider  approximating  the  modified  Bessel  Function  with  an  exponential 
function  of  the  form  f{x)  =  e"®  +13^+1  a  specified  finite  interval.  This  is  achieved  by 
fitting  a  quadratic  function,  ax^  -\-  (3x  +  '^,  to  log(/o(a^)),  where  the  coefficients  a,  [3  and  7 
are  real  constants.  Applying  the  LS  scheme,  the  coefficients  can  be  determined  using  the 
method  outlined  in  Section  1.3.  This  scheme  requires  sequential  estimation  of  the  values 
of  the  modified  Bessel  Function.  However,  this  is  not  viewed  as  a  shortcoming  because 
the  tendency  in  the  literature  is  to  construct  approximations  of  the  integral  (2)  in  terms 
of  such  functions  anyway. 


2.1.1  An  Estimator  Based  on  Integral  Representation  (6):  Ei 

In  view  of  (6),  we  can  fit  an  exponential  function  to  the  modified  Bessel  function  component 
of  the  integral.  Note  that 


J  a 


dx 


I -2a 


-f 


(2  —  4a) —  a 


Erf 


-  —  a  [  b  —  - - — 

2  V  1  -  2a 


5 


DSTO-TN-0761 


(11) 


where  S  =  e  2 


-  ^-5^" +7+ 


2-4a  and  Erf(-)  is  the  error  function 


Erf(a:)  =  /  e  ^  dt. 

Vtt  Jo 


(12) 


Hence,  applying  (11)  to  (6),  and  then  applying  the  result  to  (1),  we  arrive  at  the  HER 
estimator: 


l-e-'’'lo(62)l  +  Je-K“^+^^)/o(a6) 


^-(l-„)(a-T^)  _  (l_„)(6_^) 


+ 


1  -  2o; 

(2  —  4a) -^2  —  «  i.  \  *  2  V  1  2a 


-Erth/i-a(»-^ 


(13) 


where  the  coefficients  of  the  polynomial,  a,  fJ  and  7,  are  evaluated  numerically  using  the 
LSM  discussed  in  Section  1.3.  This  is  easily  done  in  Matlab  by  specifying  it  to  fit  a 
quadratic  expression  to  the  logarithm  of  data  points  generated  from  the  modified  Bessel 
function  of  order  zero. 


2.1.2  An  Approximation  Based  on  Integral  Representation  (5):  E2 


Next  we  consider  an  approximation  to  the  integrand  in  (5).  We  can,  as  previously,  apply 
the  LSM  to  approximate  the  modified  Bessel  Eunction  log(Io(aa:))  with  a  quadratic  curve, 
where  x  G  [a,  6] .  As  before,  note  that 


J  a 


=  e 


-7“^+7+2^  f 

J  a 


dx 


I -2a 


=-(^«)(-T^)'_,-(^«)(^-T^)' 


+ 


SP-s/n 


(2  —  4a) —  a 


Erf 


-  —  a  [  0  —  - - — 

2  V  1  -  2a 


Erfh/l-a(a-5-2_ 


(14) 
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1  2 

where  5  =  e~2“  +'i'+2-4a ,  Hence,  by  applying  (14)  to  (5),  and  by  an  application  of  the 
result  to  (1),  we  arrive  at  the  estimator: 


E2 


1 

1  +  e  “  Io{(0) 

5 

[e  (5  «)(“  lES  - 

-e  (i  “)(' 

2 

1  -  2a 

(5/3 0r 


(2  -  4a)  -  a  [ 


Erf  [  \l -  —  O'  (b - ^7— 

IV2  V  l-2a. 


Erf 


(15) 


2.2  Least  Squares  on  a  Semi-infinite  Interval: 


An  LS  estimator  is  now  proposed,  based  upon  the  original  Marcum  Q-Function  (2).  As  in 
the  previous  Subsections,  we  can  apply  a  quadratic  approximation  to  the  logarithm  of  the 
modified  Bessel  function  in  the  integrand  of  (2).  By  doing  this,  and  applying  the  result  to 
(1),  we  arrive  at  the  estimator: 


E3  — 


(5 


1  -  2a 


lit 


1  r  / 

__e  2 1  j/o(afe) 


(2  —  4a) —  a 


Erfc|,/l-a((,-^ 


(16) 


1  2  0^ 

where  5  =  6^“  and  Erfc(-)  is  the  complementary  error  function 

Erfc(x)  =  —=  /  e“*  dt. 

Jx 


(17) 


Next  we  consider  polynomial  approximations  of  the  entire  integrand,  using  the  LSM. 


2.3  Polynomial  Integrand  Approximations 


A  number  of  polynomial  approximations  are  now  examined.  In  particular,  we  apply  an 
approximation  with  a  polynomial  of  degree  5  to  the  entire  integrand  of  the  respective 
integral. 
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2.3.1  An  Estimator  Based  on  Integral  Representation  (6):  E4 

1  2 

To  begin,  we  apply  this  idea  to  the  integrand  of  (6).  In  particular,  xe~^^  /o(6x),  can  be 
approximated  by  a  fifth  order  polynomial,  f{x)  =  a^x^  +  a2X^  +  asx^  +  a4X^  +  asx  +  ag, 
where  x  G  [a,  b] .  It  can  be  shown  that 


e  2^  y  i^ocix^  +  a2X^  +  a2,x^  +  a^x^  +  a^x  +  dx  =  e  2^  [g{b)  —  g{a)]  ,  (18) 

where  g{t)  = 

Applying  (18)  to  (6),  then  (1),  results  in  the  following  estimator: 

^4  =  ^  [1  -  e-''lo(62)]  +  ie-H“^+''^)lo(a6)  -  [g{b)  -  5(a)] .  (19) 

The  coefficients  of  the  polynomial  are  estimated  by  the  LSM  in  Matlab. 

2.3.2  An  Approximation  Based  on  Integral  Representation  (8):  E, 

Using  an  argument  similar  to  that  as  in  the  construction  of  (19),  one  can  construct  the 
estimator 

^5  =  ^e-M“^+'>^)lo(a6)  +  e-5^'ff(a)  (20) 

based  upon  (8). 

2.3.3  An  Estimator  Based  on  the  Transformed  Integral  Representation: 

Eq 

Consider  the  transformation  a  =  b  =  and  v  =  ^  (see  [22]  for  details  of  this 

transformation)  of  the  integral  (8).  By  applying  the  transformation  to  (8)  and  (1),  we 
arrive  at  the  following  expression  for  the  BER: 

BER(7|<j,t)  =  ■^e“^^+^)/o(2v^)  +  [  e~^'’~^'"ho{2^/Tv)dv.  (21) 

2  Jo 

This  suggests  one  can  apply  the  LSM  on  the  interval  v  G  [0,  ?].  In  view  of  this,  the  final  LS 
estimator  we  consider  is  based  on  the  approximation  of  the  integrand,  Io{2y/rv), 

by  a  fifth  order  polynomial  using  the  LSM  on  the  interval  v  G  [0,?].  This  leads  to  the 
following  estimator: 

Eq  =  ^e"(^+‘^)/o(2v^)  +  e~'^g{(^).  (22) 

Next,  we  consider  the  Taylor  Series  Approach. 
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2.4  Taylor  Series  Approaches:  Ej 


Based  upon  Taylor  Series,  we  derive  an  alternative  estimator  of  the  BER  function  (1) 
using  some  known  series  representations  of  the  Marcum  Q-Function  and  the  zeroth  order 
modified  Bessel  Function  of  the  first  kind.  A  widely  known  double  series  expansion  ([24] 
and  [16])  of  the  standard  Marcum  Q-Function  (2)  is  given  by 


Q{a,b) 


=  e-(“'+^')  V 


/  2\/ 
I  —  I 
\2  J 

k\ 


k 

E 

3=0 


(23) 


A  well-known  series  representation  of  the  zeroth  order  modified  Bessel  function  Iq{x)  of 
the  first  kind  (see  [2]  and  [21])  is 


(24) 


Equations  (23)  and  (24)  can  be  applied  to  the  BER  equation  in  (1).  This  yields  a  series 
representation  of  the  BER  function,  involving  double  series: 


BER(7|a,6) 


Hence,  we  can  easily  construct  an  estimator  of  (1)  by  truncating  this  series  to  one  with  N 
terms.  As  N  increases  without  bound,  the  approximation  becomes  more  accurate.  Thus 
we  can  propose  the  following  estimator: 


N 

Ej  = 

k=0 


k 

E 

j=0 


V  2  / 


(26) 


The  next  section  will  examine  the  performance  of  the  six  LS  estimators  and  this  truncated 
Taylor  Series  estimator.  The  results  will  be  compared  to  that  of  ASQ. 
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3  Performance  Analysis 


We  now  examine  the  performance  of  the  seven  estimators  introduced  in  the  previous 
section.  ASQ  with  a  tolerance  of  will  be  used  as  a  benchmark  of  performance.  As 

ASQ  is  known  to  be  reliable,  it  has  been  used  by  members  of  the  Maritime  Air  Radar  Group 
in  EWRD  at  DSTO  extensively  to  estimate  intractable  integrals.  It  is  interesting  to  see 
whether  any  of  the  seven  estimators  have  the  same  level  of  accuracy  for  less  computation 
time.  Throughout  we  will  measure  7  in  decibels  (dB),  and  will  be  interested  in  values  of 
it  ranging  from  0  to  12dB,  as  considered  in  [6].  Appendix  A  contains  all  the  figures,  while 
Appendix  B  contains  tables,  from  this  analysis. 


3.1  LS  Estimators  Ei,  E2  and 

We  begin  by  analysing  the  first  three  estimators,  Ei,  E2  and  E^.  Figure  A.l  displays  a  plot 
comparing  estimates  of  Ei  and  ASQ,  while  Table  B.l  displays  corresponding  numerical 
estimates.  Included  in  the  table  are  the  absolute  and  relative  errors  of  the  estimates,  when 
compared  to  ASQ.  The  estimates  are  obtained  by  fitting  five  points  to  the  modified  Bessel 
Function.  These  results  show  that  Ei  performs  well  on  the  region  7  <  1,  with  increasing 
accuracy  for  small  SNR  values. 

Next  we  examine  E2.  Figure  A. 2  provides  a  plot  of  E2  in  comparison  to  ASQ.  Table 
B.2  provides  the  approximations  of  BFR  using  E2.  Observe  that  E2  outperforms  Ei  for 
extremely  small  values  of  7,  more  specifically,  for  7  <  3.  However,  the  accuracy  worsens 
dramatically  for  larger  values  of  7.  Beyond  7  >  8,  F2  becomes  relatively  inaccurate,  and 
thus  results  are  not  included  in  Table  B.2.  The  approximation  in  E2  performs  well  only 
on  the  interval  where  7  is  extremely  small. 

The  performance  of  F3  can  be  viewed  in  Figure  A. 3,  with  numerical  results  in  Table  B.3.  It 
is  interesting  to  note  that  the  estimators  Ei  and  E2  are  generally  superior  to  estimator  F3 
for  small  values  of  7,  but  their  performance  deteriorates  quickly.  In  contrast,  E^  provides 
better  approximations  for  larger  values  of  7,  and  the  accuracy  improves  as  7  increases. 
However,  these  three  estimators  have  relatively  large  errors  and  are  not  entirely  consistent 
with  ASQ.  The  accuracy  can  be  slightly  improved  by  increasing  the  number  of  points 
being  fitted  to  the  modified  Bessel  Function. 

As  a  final  comparison  of  these  three  estimators,  Figure  A. 4  shows  all  on  the  same  plot, 
with  ASQ. 


3.2  LS  Estimators  E4,  E^  and  Eq 

We  now  examine  the  estimators  E4,  E^  and  Eq.  Figure  A. 5  provides  a  plot  comparing 
E4  and  ASQ.  In  Table  B.4,  estimates  of  BFR  using  E4  are  compared  with  those  obtained 
via  ASQ.  The  actual  error  and  relative  error  are  given,  as  previously.  One  hundred  points 


10 


DSTO-TN-0761 


equally  spaced  on  the  interval  [a,  b]  are  used  for  the  LS  computation.  As  can  be  observed, 
this  estimator  gives  better  results  than  the  estimators  Ei,  E2  and  E^  for  small  values  of  7, 
with  increasing  accuracy  for  smaller  SNR  values.  However,  the  estimator’s  performance 
deteriorates  rapidly  for  larger  values  of  7.  To  improve  the  accuracy,  one  can  increase  the 
order  of  the  polynomial  and  the  number  of  points  being  fit  to  the  modified  Bessel  Function. 

Figure  A. 6  is  a  plot  of  the  performance  of  estimator  E^.  As  can  be  seen,  there  is  a  uniform 
improvement  in  performance.  Table  B.5  confirms  the  improvement  in  terms  of  errors.  We 
observe  that  E^  is  generally  superior  to  E4.  The  results  are  accurate  for  a  larger  range  of 
7  values.  Although  the  accuracy  of  the  approximation  decreases  as  7  increases,  the  rate 
of  deterioration  is  relatively  slower  than  all  estimators  considered  thus  far. 

We  now  consider  the  final  LS  estimator,  Eq.  Figure  A. 7  shows  its  performance  relative  to 
ASQ.  This  estimator  proved  to  be  the  most  accurate  and  efficient  LS  estimator.  Included  in 
Figure  A. 7  are  estimates  derived  from  the  results  in  [6],  based  upon  (3).  Table  B.6  contains 
results  using  Eq,  (3)  and  ASQ.  Estimator  Eq  used  50  points  for  fitting  the  modified  Bessel 
function.  Relative  errors  are  given,  with  ei  that  between  Eq  and  ASQ,  and  62  the  relative 
error  given  in  [6].  As  shown  in  Table  B.6,  Eq  outperforms  (3)  and  other  LS  estimators 
discussed  so  far.  It  is  worth  noting  that  Eq  provides  extremely  accurate  results  for  the 
entire  range  of  7  values  of  interest.  Notice  that  the  actual  error  is  relatively  consistent. 
As  for  7  larger  than  12dB,  the  proposed  approximation  (3)  in  [6]  remains  a  better  choice. 
Eq  may  be  of  practical  interest  when  7  is  less  than  12dB,  which  as  pointed  out  in  [6],  is 
the  region  of  interest. 

Figure  A. 8  shows  the  three  estimators  considered  in  this  subsection  in  the  same  plot,  with 
ASQ. 

In  addition,  we  experiment  with  the  degree  of  the  approximating  polynomial,  n,  in  Eq.  For 
comparison,  we  include  the  results  of  Eq  using  n  =  6  with  (3)  and  (4)  in  [6]  on  the  interval 
7  between  OdB  and  IddR  in  Table  B.7  .  The  corresponding  relative  errors  are  shown.  Note 
that  Eq  using  n  =  6  provides  better  results,  with  accuracy  as  high  as  2  x  10~^^  for  some 
values  of  7.  The  approximations  are  of  extremely  high  accuracy  for  the  region  of  interest. 
Figure  A. 9  provides  a  corresponding  comparison  plot.  Figure  A.  10  shows  an  enlarged 
image  of  Figure  A. 9  around  7  =  13.  The  efficiency  performance  of  Eq  is  analyzed  in  the 
subsequent  subsection.  It  is  important  to  note  that  Eq  is  capable  of  improved  accuracy 
results  by  increasing  the  value  of  n,  namely  the  degree  of  the  LS  polynomial.  One  can 
vary  the  value  of  n  to  obtain  the  desired  accuracy. 

This  completes  our  examination  of  LS  estimators.  We  now  investigate  the  Taylor  Series 
estimator,  Ej. 


3.3  Taylor  Series  Estimator  Ej 

The  final  estimator  we  consider  is  Ej,  which  is  based  upon  the  truncated  series  in  (26). 
The  partial  sum  uses  80  terms.  Figure  A. 11  contains  estimates  of  Ej.  It  is  worth  noting 
that  the  estimates  given  by  (26)  are  extremely  accurate,  especially  for  small  values  of  7.  It 
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clearly  outperforms  the  LS  estimators  given  in  the  previous  subsections.  We  do  not  include 
a  comparison  plot  between  Eq  and  E^  because  the  errors  in  Table  B.8,  when  compared  to 
the  results  in  Table  B.6,  indicate  an  enormous  improvement.  Note  that,  for  larger  values 
of  7,  the  relative  error  increases,  but  the  accuracy  can  be  further  improved  by  increasing 
N,  the  number  of  terms  in  the  partial  sum. 

We  now  consider  computational  times  associated  with  these  estimators.  We  will,  in  par¬ 
ticular,  be  interested  in  how  well  Eq  and  Ej  perform  relative  to  ASQ. 


3.4  Time  Performance  Analysis 

The  previous  subsections  identified  estimators  Eq  and  Ej  as  having  the  most  consistent 
performance  with  ASQ,  and  in  particular,  Ej  is  the  more  accurate  of  the  two.  It  will  thus 
be  useful  to  investigate  processing  speeds  associated  with  these  estimators.  Specifically, 
we  will  be  interested  in  knowing  whether  either  one  has  comparable  performance  to  ASQ 
in  shorter  processing  times. 

The  time  performance  of  Ej  is  compared  with  Eq  and  ASQ  in  Figures  A. 12  and  A. 13,  as 
well  as  in  Table  B.9.  Computation  time  is  measured  in  seconds,  while  accuracy  is  specified 
in  negative  powers  of  10.  Note  that  in  view  of  Table  B.9,  Eq  is  the  most  efficient  estimator, 
while  Ej  is  the  second  and  ASQ  the  least.  Moreover,  Eq  and  E’r  are  very  consistent.  At 
7  =  OdB,  Eq  obtained  an  exact  result  with  a  polynomial  of  degree  n  =  8,  while  E^  with  a 
Taylor  Series  with  partial  sum  N  =  12.  In  contrast,  for  ASQ,  time  complexity  increases  as 
the  accuracy  increases.  A  plot  of  the  results  in  Table  B.9  is  given  in  Figure  A.  12.  Figure 
A. 13  gives  an  enlarged  view  of  Figure  A. 12  at  an  accuracy  of  10“^^.  Notice  that  Eq  slightly 
outperforms  £'7. 
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4  Conclusion 


This  report  investigated  numerical  estimates  of  the  BER  function,  in  particular,  via  LS  and 
Taylor  Series  approximations  on  the  Marcum  Q-Function.  Using  a  number  of  finite  integral 
representations  of  the  Marcum’s  Q-Function,  several  LS  estimators  were  derived.  Three 
were  based  on  approximating  the  modified  Bessel  function  by  an  exponential  function  with 
quadratic  argument  {Ei,E2,  E2).  Another  three  were  based  upon  a  fifth  order  polynomial 
approximation  of  the  entire  integrand  {E/^,E^,Eq).  An  estimator,  based  upon  Taylor 
Series  was  also  introduced  (Ej). 

Their  accuracy  and  efficiency  performances  were  analysed.  All  LS  estimators  perform 
well  on  certain  local  regions.  Results  indicated  that  the  estimators,  E4,  E5  and  Eq,  were 
generally  superior  to  Ei,  E2  and  F3.  The  optimal  LS  estimator  identified  in  this  report  is 
Eq,  providing  the  highest  accuracy  on  the  entire  region  of  interest.  It  is  the  most  efficient 
estimator,  with  higher  efficiency  than  the  Taylor  Series  estimator  Ej. 

From  the  accuracy  perspective,  Ey  outperforms  the  LS  estimators  derived  in  this  report. 
High  accuracy  results  were  obtained  at  fast  speed.  It  is  worth  noting  that  the  series 
approximation  requires  much  less  computation  time  than  the  ASQ  approach  to  achieve 
the  same  accuracy. 

The  results  presented  here  also  demonstrated  that  the  estimates  (3)  and  (4)  from  [6]  can 
be  improved  significantly  within  the  region  of  interest  (0  to  12dB),  by  using  Eq,  but  for 
estimates  greater  than  12dB,  the  estimate  (4)  is  suitable. 

It  is  worth  noting  that  the  Least  Squares  Method  is  already  available  in  the  computer 
language  Matlab,  and  so  the  general  methodology  used  here  may  be  applied  to  estimate 
other  intractable  integrals  of  interest. 


5  Acknowledgements 


The  authors  would  like  to  thank  the  sponsor.  Commander  Surveillance  &  Response  Group, 
RAAF  Williamtown,  for  supporting  this  research  through  Task  AIR  04/206.  The  Co¬ 
author  Sharon  Lee  expresses  thanks  to  DSTO  for  providing  the  opportunity  to  be  involved 
in  its  Summer  Vacation  Scholarship  Program  06/07.  Thanks  are  due  to  Dr  Paul  Berry  for 
vetting  the  report. 


13 


DSTO-TN-0761 


6  References 

1.  Breed,  G.,  Bit  Error  Rates:  Fundamental  Concepts  and  Measurement  Issues,  High  Freq. 
Elec.  1,  46-48,  2003. 

2.  Bowman,  F.,  Introduction  to  Bessel  Funetions.  (Dover,  New  York,  1953). 

3.  Burden,  R.  L.  and  Faires,  J.  D.,  Numerical  Analysis.  (Brooks  Cole,  US,  2004). 

4.  Chiani,  M.,  Integral  representation  and  bounds  for  Marcum  Q-Function,  IFF  Flee.  Lett. 
35,  445-446,  1999. 

5.  Ferrari,  G.  and  Corazza,  G.E.,  New  bounds  on  the  Marcum  Q-Function,  IFFF  Trans. 
Inf.  Theory  48,  3003-3008,  2002. 

6.  Ferrari,  G.  and  Corazza,  G.E.,  Tight  bounds  and  accurate  approximations  for  DQPSK 
transmission  bit  error  rate,  IFF  Flee.  Lett.  40,  1284-1285,  2004. 

7.  Kam,  P.  Y.  and  Li,  R.,  A  New  Geometric  View  of  the  First-Order  Marcum  Q-Function 
and  Some  Simple  Tight  Erfe-Bounds,  IFFF  63rd  Veh.  Tech.  Conf.  5,  2553-2557,  2006. 

8.  Kam,  P.  Y.  and  Li,  R.,  Simple  Tight  Exponential  Bounds  on  the  First-Order  Marcum 
Q-Function  via  the  Geometric  Approach,  IFFF  Int.  Sym.  Trans.  Info.  Theory  7,  1085- 
1089,  2006. 

9.  Kam,  P.  Y.  and  Li,  R.,  Computing  and  Bounding  the  Generalized  Marcum  Q-Function 
via  a  Geometric  Approach,  IFFF  Int.  Sym.  Trans.  Info.  Theory  7,  1090-1094,  2006. 

10.  Lyness,  J.  N.  and  Kaganove,  J.  J.,  Comments  on  the  Nature  of  Automatic  Quadrature 
Routines.  ACM  Trans.  Math.  Soft.  2  65-81,  1976. 

11.  Marcum,  J.  L,  Tables  of  Q  functions.  Rand  Corp.,  Santa  Monica,  CA.  U.S.A.F.  Project 
RAND  Research  Memorandum  M-339,  1950. 

12.  Marcum,  J.  L.,  A  statistical  theory  of  target  detection  by  pulsed  radar.  IRF  Trans. 
Info.  Theory,  IT-6,  59-144,  1960. 

13.  Marcum,  J.  1.  and  Swerling,  P.,  Studies  of  target  detection  by  pulsed  radar.  IFFF 
Trans.  Info.  Theory,  IT-6,  1960. 

14.  Proakis,  J.  G.,  Digital  Communications.  (McGraw-Hill,  New  York,  2001). 

15.  Schwartz,  M.,  Bennett,  W.  R.  and  Stein,  S.  (1996),  Communication  Systems  and 
Techniques.  (McGraw-Hill,  New  York,  1996). 

16.  Shnidman,  D.  A.,  The  Calculation  of  the  Probability  of  Detection  and  the  Generalized 
Marcum  Q-Function.  IFFF  Trans.  Theory  35,  389-400,  1989. 

17.  Simon,  M.  K.,  A  New  Twist  on  the  Marcum  Q-Function  and  its  Application,  IFFF 
Comms.  Lett.,  2,  39-41,  1998. 

18.  Simon,  M.K.  and  Alouini,  M.  S.,  A  Unified  Approach  to  the  Probability  of  Error  for 
Noncoherent  and  Differentially  Coherent  Modulations  Over  Generalized  Fading  Chan¬ 
nels,  IFFF  Trans.  Comm.  46,  1625-1638,  1998. 


14 


DSTO-TN-0761 


19.  Simon,  M.K.  and  Alouini,  M.  S.,  Exponential-type  bounds  on  the  generalized  Marcum 
Q-Function  with  application  to  error  probability  analysis  over  fading  channels,  IEEE 
Trans.  Comms.,  3,  359-366,  2000. 

20.  Simon,  M.  K.  and  Alouini,  M.  S.,  Some  New  Results  for  Integrals  involving  the  Gen¬ 
eralised  Marcum  Q-  Function  and  their  Application  to  Performance  Evaluation  Over 
fading  Channels.  IEEE  Trans.  Wireless  Comms.,  2,  611-615,  2003. 

21.  Tsypkin,  A.  G.  and  Tsypkin,  G.  G.,  Mathematical  Formulas:  Algebra,  Geometry  and 
Mathematical  Analysis.  (Mir  Publishers,  Moscow,  1988). 

22.  Weinberg,  G.V.,  Stochastic  Representations  of  the  Marcum  Q-Function  and  Associ¬ 
ated  Radar  Detection  Probabilities.  DSTO-RR-0304,  2005. 

23.  Weinberg,  G.  V.  and  Panton,  L.,  Numerical  Estimation  of  Marcum’s  Q-Function  using 
Monte  Carlo  Approximation  Schemes.  DSTO-RR-0311,  2006. 

24.  Weinberg,  G.  V.,  A  Poisson  Representation,  and  Monte  Carlo  Estimation,  of  the 
Generalized  Marcum  Q-Function,  IEEE  Trans.  Aero.  Elec.  Sys.  42,  1520-1531,  2006. 


15 


DSTO-TN-0761 


16 


DSTO-TN-0761 


Appendix  A:  Comparsion  Plots  of  the  Estimators 


Plot  of  the  BER 


Figure  A.l:  A  plot  comparing  Ei  and  ASQ  for  the  region  7  =  OdB  to  ^  =  12dB.  The 
BER  is  given  in  a  logarithmic  scale.  Observe  that  estimator  Ei  performs  well  for  small 
values  of'j,  specifically  where  7  <  QdB. 
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Plot  of  the  BER 


Figure  A. 2:  Performance  of  E2  in  comparision  with  ASQ  on  the  interval  7  =  OdB  to  ^  = 
8dB.  The  BER  is  given  in  a  logarithmic  scale.  Estimator  E2  provides  good  approximation 
for  7  <  bdB,  with  increasing  accuracy  for  smaller  SNR  values. 


Plot  of  the  BER 


Figure  A. 3:  A  plot  displaying  the  performance  of  E3  in  comparision  to  ASQ  for  7  =  OdB 
to  ^  =  12dB.  The  BER  are  given  in  a  logarithmic  scale.  Estimator  E^  performs  well  for 
larger  values  of  'j,  in  particular,  7  >  lOdB.  Observe  that  accuracy  increases  as  the  SNR 
increases. 
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Plot  of  the  BER 


Figure  A. 4:  A  Comparison  of  the  three  estimators  Ei,  E2  and  E3,  with  ASQ. 


Plot  of  BER  in  log  scale 


Figure  A. 5:  Performance  0/ in  comparision  with  ASQ  on  the  interval  7  =  ddB  to 
7  =  lOdS.  The  BER  is  given  in  a  logarithmic  scale.  Notice  that  for  7  beyond  8dB, 
Ei  begins  to  deviate  from  ASQ.  Estimator  E^  provides  good  approximation  for  the  region 
7  <  8dB,  with  increasing  accuracy  for  smaller  SNR  values. 
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Plot  of  BER  in  log  scale 


Figure  A. 6:  A  plot  comparing  E5  and  ASQ  for  the  region  7  =  OdB  to  'y  =  12dB.  Observe 
that  E5  covers  ASQ.  For  the  range  of'j  values  of  interest,  approximation  E^  provides  good 
results.  High  accuracy  approximations  are  obtained  for  7  <  9dB. 


Plot  of  BER  in  log  scale 


Figure  A.  7:  Performance  of  Eq  in  comparision  with  ASQ  on  the  interval  7  =  9dB  to 
7  =  l2dB.  The  BER  is  given  in  a  logarithmic  scale.  Approximation  Eq  performs  extremely 
well  on  the  region  of  interest,  providing  high  accuracy  results  for  7  <  12dB. 


20 


DSTO-TN-0761 


Plot  of  BER  in  log  scale 


Figure  A. 8:  A  Comparison  of  the  three  estimators  E^,  and  Eq,  with  ASQ. 


Plot  of  BER  in  log  scale 


Figure  A. 9:  A  comparison  plot  of  Eq  using  n  =  6,  (3),  (4)  and  ASQ  for  7  =  QdB  to 
7  =  lAdB.  The  BER  is  given  in  a  logarithmic  scale.  Estimator  Eq  with  n  =  6  performs 
extremely  well  on  this  region,  providing  extremely  high  accuracy  results.  Observe  that  Eq 
with  n  =  6  is  almost  exactly  on  ASQ. 
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Plot  of  BER  in  log  scale 


Figure  A. 10:  An  enlarged  view  of  Figure  A. 9  around  7  =  13.  Note  that  BER  using  (4) 
falls  outside  the  region  of  the  Figure.  Observe  how  closely  Eq  approximate  ASQ. 


Plot  of  BER  in  log  scale 


Figure  A. 11:  A  plot  of  Ej  and  ASQ  for  7  =  OdB  to  7  =  12dB.  The  BER  is  given  in  a 
logarithmic  scale.  Extremely  high  accuracy  results  is  obtained  using  the  series  in  (26)  for 
the  entire  region  of  interest. 
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Performance  measure  (time)  of  ASQ,  E6  and  E7  at  OdB 
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Figure  A.  12:  Time  performance  plot  of  Eq,  Ej  and  ASQ.  Time  is  in  seconds  and  accuracy 
is  given  in  a  logarithmic  scale,  base  1 0.  Observe  that  both  Eq  and  E-j  are  extremely  efficient 
in  comparision  with  ASQ.  Both  Eq  and  E'y  are  very  consistent,  while  ASQ  depends  heavily 
on  the  degree  of  accuracy. 


X  10^^  Performance  measure  (time)  of  ASQ,  E6  and  E7  at  OdB 
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Figure  A.  13:  An  enlarged  view  of  Figure  A.  12  at  10  suggesting  Eq  is  more  efficient 
than  E-;.  Note  that  ASQ  falls  outside  the  region  of  the  Figure. 
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Appendix  B:  Tables  of  Numerical  Results 


l{dB) 

ASQ 

£1 

ei 

€2 

0 

1.63907530399585e-001 

1.64066250992386e-001 

1.5872e-004 

9.6835e-002% 

1 

1.30332229277756e-001 

1.30554573532828e-001 

2.2234e-004 

1.7060e-001% 

2 

9.93324252133400e-002 

9.96169239202047e-002 

2.8450e-004 

2.8641e-001% 

3 

7.18727612659082e-002 

7.22160657416420e-002 

3.4330e-004 

4.7766e-001% 

4 

4.87488622380308e-002 

4.91535216927732e-002 

4.0466e-004 

8.3009e-001% 

5 

3.04943244285625e-002 

3.09694472029306e-002 

4.7512e-004 

1.5581e+000% 

6 

1.72359006046928e-002 

1. 7791392359341 le-002 

5.5549e-004 

3.2229e+000% 

7 

8.58004434130124e-003 

9.22092893335622e-003 

6.4088e-004 

7.4695e+000% 

8 

3.64294312896472e-003 

4.36665001049495e-003 

7.2371e-004 

1.9866e+001% 

9 

1.26710356868390e-003 

2.06276280461294e-003 

7.9566e-004 

6.2794e+001% 

10 

1.19230674177357e-003 

3.43184596033452e-004 

8.4912e-004 

2.4742e+002% 

11 

9.46488998159275e-004 

9.46488998159275e-004 

8.7859e-004 

1.2940e+003% 

12 

9.05258912217385e-006 

8.90723374721225e-004 

8.8167e-004 

9.7394e+003% 

Table  B.l:  A  comparision  of  Ei  and  ASQ  approximations,  with  a  tolerance  o/10~^®.  Five 
equally  spaced  points  are  used  for  fitting  the  modified  Bessel  function.  ei  represents  the 
absolute  error  between  Ei  and  ASQ,  and  €2  is  the  relative  error 


l{dB) 

ASQ 

£2 

ei 

€2 

0 

1.63907530399585e-001 

1.63896182918470e-001 

1.1347e-005 

6.9231e-003% 

1 

1.30332229277756e-001 

1.30283201772424e-001 

4.9028e-005 

3.7617e-002% 

2 

9.93324252133400e-002 

9.91910467676262e-002 

1.4138e-004 

1.4233e-001% 

3 

7.18727612659082e-002 

7.15574732746942e-002 

3.1529e-004 

4.3868e-001% 

4 

4.87488622380308e-002 

4.81818860883209e-002 

5.6698e-004 

1.1631e+000% 

5 

3.04943244285625e-002 

2.96511574388189e-002 

8.4317e-004 

2.7650e+000% 

6 

1.72359006046928e-002 

1.61643524537770e-002 

1.0715e-003 

6.2170e+000% 

7 

8.58004434130124e-003 

7.35773400350712e-003 

1.2223e-003 

1.4246e+001% 

8 

3.64294312896472e-003 

2.31931668012224e-003 

1.3236e-003 

3.6334e+001% 

Table  B.2:  Approximations  of  BER{'y\a, 
to  'y  =  8dB.  Five  equally  spaced  points 
ei  =  \  ASQ  —  £’2!  and  €2  =  100  x  ^4^. 


b)  based  on  estimator  £2  on  the  interval  7  =  OdB 
are  used  for  fitting  the  modified  Bessel  function. 
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l{dB) 

ASQ 

£3 

ei 

€2 

0 

1.63907530399585e-001 

1.71550949650647e-001 

7.6434e-003 

4.6633e+000% 

1 

1.30332229277756e-001 

1.37826804088733e-001 

7.4946e-003 

5.7504e+000% 

2 

9.93324252133400e-002 

1.05719011266452e-001 

6.3866e-003 

6.4295e+000% 

3 

7.18727612659082e-002 

7.64823916401551e-002 

4.6096e-003 

6.4136e+000% 

4 

4.87488622380308e-002 

5.15201561338994e-002 

2.7713e-003 

5.6848e+000% 

5 

3.04943244285625e-002 

3.18780295476395e-002 

1.3837e-003 

4.5376e+000% 

6 

1.72359006046928e-002 

1. 781651 19904009e-002 

5.8061e-004 

3.3686e+000% 

7 

8.58004434130124e-003 

8.78891826423036e-003 

2.0887e-004 

2.4344e+000% 

8 

3.64294312896472e-003 

3.70801023832913e-003 

6.5067e-005 

1.7861e+000% 

9 

1.26710356868390e-003 

1.28431703460022e-003 

1.7213e-005 

1.3585e+000% 

10 

1.19230674177357e-003 

3.46862265107066e-004 

3.6777e-006 

1.0716e+000% 

11 

9.46488998159275e-004 

6.84867934099513e-005 

5.9167e-007 

8.7144e-001% 

12 

9.05258912217385e-006 

9.11841944014908e-006 

6.5830e-008 

7.2720e-001% 

Table  B.3:  A  comparision  of  with  ASQ  approximations.  Five  equally  spaced  points  are 
used  for  fitting  the  modified  Bessel  function.  ei  represents  the  relative  error  of  and  €2 
is  the  relative  error. 


l{dB) 

ASQ 

£4 

ei 

€2 

0 

1.63907530399585e-001 

1.63907488870700e-001 

4.1529e-008 

2.5337e-005% 

1 

1.30332229277756e-001 

1.30332106347130e-001 

1.2293e-007 

9.4321e-005% 

2 

9.93324252133400e-002 

9.93320691796630e-002 

3.5603e-007 

3.5843e-004% 

3 

7.18727612659082e-002 

7.18718252819206e-002 

9.3598e-007 

1.3023e-003% 

4 

4.87488622380308e-002 

4.87466026625585e-002 

2.2596e-006 

4.6351e-003% 

5 

3.04943244285625e-002 

3.04892302562418e-002 

5.0942e-006 

1.6705e-002% 

6 

1.72359006046928e-002 

1.72251622791136e-002 

1.0738e-005 

6.2302e-002% 

7 

8.58004434130124e-003 

8.55909104642799e-003 

2.0953e-005 

2.4421e-001% 

8 

3.64294312896472e-003 

3.60577177556770e-003 

3.7171e-005 

1.0204e-h000% 

9 

1.26710356868390e-003 

1.20888066573677e-003 

5.8223e-005 

4.5950e-h000% 

10 

3.43184596033453e-004 

2.67003180762893e-004 

7.6181e-005 

2.2198e+001% 

11 

6.78951281813439e-005 

-3.6969541 1045479e-006 

7.1592e-005 

1.0545e-h002% 

12 

9.05258912217385e-006 

-2.75618525530374e-006 

1.1809e-005 

1.3045e-h002% 

Table  B.4:  Estimates  of  BE R{'y\a,b)  using  £'4  on  the  interval  7 
100  equally  spaced  points  are  used  for  fitting  the  intergand  in  (6). 


€2  =  100  X 


£1 

ASQ- 


=  OdB  to  7  =^2dB. 
€i  =  \ASQ  —  Ef\  and 
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l{dB) 

ASQ 

A5 

ei 

€2 

0 

1.63907530399585e-001 

1.63907546749063e-001 

1.5105e-008 

9.2154e-006% 

1 

1.30332229277756e-001 

1.30332244210175e-001 

1.3768e-008 

1.0564e-005% 

2 

9.93324252133400e-002 

9.93324155663120e-002 

1.0629e-008 

1.0700e-005% 

3 

7.18727612659082e-002 

7.18726865834308e-002 

7.5848e-008 

1.0553e-004% 

4 

4.87488622380308e-002 

4.87486871417470e-002 

1.7599e-007 

3.6102e-004% 

5 

3.04943244285625e-002 

3.04940693189189e-002 

2.5575e-007 

8.3867e-004% 

6 

1.72359006046928e-002 

1.72356861246986e-002 

2.1509e-007 

1.2479e-003% 

7 

8.58004434130124e-003 

8.58005712697456e-003 

7.3343e-009 

8.5480e-005% 

8 

3.64294312896472e-003 

3.64328283674178e-003 

3.3717e-007 

9.2554e-003% 

9 

1.26710356868390e-003 

1.26763628874464e-003 

5.3069e-007 

4.1882e-002% 

10 

3.43184596033453e-004 

3.43637039076888e-004 

4.4649e-007 

1.3010e-001% 

11 

6.78951281813439e-005 

6.81247127723810e-005 

2.2577e-007 

3.3250e-001% 

12 

9.05258912217385e-006 

9.12061302541548e-006 

6.7384e-008 

7.4431e-001% 

Table  B.5:  A  comparision  of  the  approximation  of  BER{'y\a,b)  based  on  equation  (20) 
with  that  obtained  via  ASQ.  50  equally  spaced  points  are  used  for  fitting  the  intergrand  in 
(8).  ei  =  \ASQ  —  E^\  and  62  =  100  x 


AdB) 

ASQ 

Eq 

ei 

BER  in  [6] 

€2 

0 

1.63907530399585e-001 

1.63907530400071e-001 

2.9656e-010% 

1.73998678128697e-001 

6.1566e+000% 

1 

1.30332229277756e-001 

1.30332229279850e-001 

1.6061e-009% 

1.36604369075718e-001 

4.8124e+000% 

2 

9.93324252133400e-002 

9.93324252171851e-002 

3.8709e-009% 

1.02741930798427e-001 

3.4324e+000% 

3 

7.18727612659082e-002 

7.18727612659082e-002 

1.2320e-008% 

7.35176792703302e-002 

2.2887e+000% 

4 

4.87488622380308e-002 

4.87488621670844e-002 

1.4553e-007% 

4.94684302616405e-002 

1.4761e+000% 

5 

3.04943244285625e-002 

3.04943243033633e-002 

4.1057e-007% 

3.07841984642038e-002 

9.5058e-001% 

6 

1.72359006046928e-002 

1.72359008025161e-002 

1.1477e-006% 

1.73427752950751e-002 

6.2007e-001% 

7 

8.58004434130124e-003 

8.58004526916485e-003 

1.0814e-005% 

8.61510874905393e-003 

4.0867e-001% 

8 

3.64294312896472e-003 

3.64294348167390e-003 

9.6820e-006% 

3.65278932995171e-003 

2.7028e-001% 

9 

1.26710356868390e-003 

1.26710160675231e-003 

1.5484e-004% 

1.26936632743277e-003 

1.7858e-001% 

10 

3.43184596033453e-004 

3.43182757950219e-004 

5.3560e-004% 

3.43588251021892e-004 

1.1762e-001% 

11 

6.78951281813439e-005 

6.78960539104487e-005 

1.3635e-003% 

6.79475127646939e-005 

7.7155e-002% 

12 

9.05258912217385e-006 

9.05419050617861e-006 

1.7690e-002% 

9.05715038668292e-006 

5.0386e-002% 

Table  B.6:  Performance  of  Eq,  (3)  and  ASQ.  ei  represents  the  relative  error  of  Eq  and  €2 
is  the  relative  error  given  in  [5],  Observe  that  E^  outperforms  (3). 
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'y(dB) 

ASQ 

Eq  with  n  =  6 

€1 

BER  using  (3) 

£2 

BER  using  (4) 

€3 

0 

1.639075303e-001 

1.639075303e-001 

1.6934e-014% 

1.739986781e-001 

6.1566e+000% 

2.254210348e-001 

3.7529e+001% 

1 

1.303322292e-001 

1.303322292e-001 

1.4907e-013% 

1.366043690e-001 

4.8124e+000% 

1.726326057e-001 

3.2456e+001% 

2 

9.933242521e-002 

9.933242521e-002 

8.6621e-013% 

1.027419307e-001 

3.4324e+000% 

1.271145833e-001 

2.7969e+001% 

3 

7.187276126e-002 

7.187276126e-002 

1.1334e-011% 

7.351767927e-002 

2.2887e+000% 

8.908322004e-002 

2.3946e+001% 

4 

4.874886223e-002 

4.874886223e-002 

1.6426e-011% 

4.946843026e-002 

1.4761e+000% 

5.866306626e-002 

2.0337e+001% 

5 

3.049432442e-002 

3.049432442e-002 

4.7817e-010% 

3.078419846e-002 

9.5058e-001% 

3.571928426e-002 

1.7134e+001% 

6 

1.723590060e-002 

1.723590060e-002 

2.1835e-009% 

1.734277529e-002 

6.2007e-001% 

1.970599722e-002 

1.4331e+001% 

7 

8.580044341e-003 

8.580044339e-003 

1.5631e-008% 

8.615108749e-003 

4.0867e-001% 

9.601935885e-003 

1.1910e+001% 

8 

3.642943128e-003 

3.642943125e-003 

9.1022e-008% 

3.652789329e-003 

2.7028e-001% 

4.001461897e-003 

9.8415e+000% 

9 

1.267103568e-003 

1.267103576e-003 

5.8926e-007% 

1.269366327e-003 

1.7858e-001% 

3.431846049e-004 

8.0898e+000% 

10 

3.43184596e-004 

3.431846049e-004 

2.5988e-006% 

3.435882510e-004 

1.1762e-001% 

7.155565315e-005 

6.6184e+000% 

11 

6.78951116e-005 

6.789511160e-005 

2.4416e-005% 

6.794751276e-005 

7.7155e-002% 

7.155565315e-005 

5.3914e+000% 

12 

9.05258912e-006 

9.052581472e-006 

8.4507e-005% 

9.057150386e-006 

5.0386e-002% 

9.448643705e-006 

4.3750e+000% 

13 

7.34977835e-007 

7.349868908e-007 

1.2321e-003% 

7.352186001e-007 

3.2758e-002% 

7.609826574e-007 

3.5382e+000% 

14 

3.19776721e-008 

3.198172372e-008 

1.2670e-002% 

3.198445315e-008 

2.1205e-002% 

3.288995544e-008 

2.8529e+000% 

Table  B.7:  A  comparision  of  Eq  with  n  =  6,  approximations  using  [6]  and  ASQ  on  the 
interval  7  between  OdB  and  lAdB.ei  represents  the  relative  error  of  Eq  with  n  =  6,  62 
indicates  the  relative  error  of  (3)and  €3  is  the  relative  error  of  (4)-  Notice  the  high  degree 
of  accuracy  of  Eq  for  small  values  0/7. 


l{dB) 

ASQ 

E^ 

ei 

€2 

0 

1.63907530399585e-001 

1.63907530399585e-001 

O.OOOOe+000 

O.OOOOe+000% 

1 

1.30332229277756e-001 

1.30332229277756e-001 

O.OOOOe+000 

O.OOOOe+000% 

2 

9.93324252133400e-002 

9.93324252133400e-002 

O.OOOOe+000 

O.OOOOe+000% 

3 

7.18727612659082e-002 

7.18727612659082e-002 

4.1633e-017 

5.7926e-016% 

4 

4.87488622380308e-002 

4.87488622380307e-002 

5.5511e-017 

1.1387e-015% 

5 

3.04943244285625e-002 

3.04943244285625e-002 

1.0408e-017 

3.4132e-016% 

6 

1.72359006046928e-002 

1.72359006046928e-002 

1.7347e-017 

1.0065e-015% 

7 

8.58004434130124e-003 

8.58004434130123e-003 

5.2042e-018 

6.0654e-016% 

8 

3.64294312896472e-003 

3.64294312896472e-003 

4.3368e-019 

1.1905e-016% 

9 

1.26710356868390e-003 

1.26710356868390e-003 

6.5052e-019 

5.1339e-016% 

10 

3.43184596033453e-004 

3.43184596033452e-004 

3.7947e-019 

1.1057e-015% 

11 

6.78951281813439e-005 

6.78951281813436e-005 

2.8460e-019 

4.1918e-015% 

12 

9.05258912217385e-006 

9.05258912217359e-006 

2.5919e-019 

2.8632e-014% 

Table  B.8:  An  approximation  of  BE R{^\a,b) 
eguation  (26). ei  =  \ASQ  —  Er\  and  62  =  100  x 


based  on  a  partial  sum  of  80  terms  using 

ei 

ASQ- 
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Accuracy 

ASQ 

Eq 

Er 

i(r^ 

3.7900e-002 

3.7935e-003 

3.6175e-003 

10-9 

4.6300e-002 

3.7935e-003 

4.3964e-003 

10-10 

6.1700e-002 

3.7935e-003 

5.4035e-003 

10-11 

8.6800e-002 

6.5642e-003 

5.4035e-003 

10-12 

1.3660e-001 

6.5642e-003 

6.6447e-003 

10-13 

2.1310e-001 

6.5642e-003 

6.6447e-003 

10-14 

4.3790e-001 

5.7586e-003 

8.7316e-003 

10-15 

4.8680e-001 

5.7586e-003 

1.0030e-002 

10-16 

8.6340e-001 

5.7586e-003 

1.0030e-002 

10-17 

1.2287e+000 

5.7586e-003 

1.2531e-002 

10-18 

1.9540e+000 

5.8133e-003 

1.2531e-002 

Table  B.9:  Time  performance  of  Eq,  Ej  and  ASQ.  Time  is  given  in  seconds.  Accuracy  is 
given  by  where  m  G  {8,  9, . . . ,  18}.  Observe  that  Eq  is  the  most  efficient  estimator, 

with  Ej  following  closely  behind.  ASQ  is  the  least  efficient,  and  time  increases  as  the 
accuracy  increases.  Both  Eq  and  Ej  are  very  consistent. 
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